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In a previous work, we introduced a connputational nnodel of area 3b which is built upon 
the neural field theory and receives input fronn a simplified nnodel of the index distal finger 
pad populated by a random set of touch receptors (Merkell cells). This model has been 
shown to be able to self-organize following the random stimulation of the finger pad model 
and to cope, to some extent, with cortical or skin lesions. The main hypothesis of the 
model is that learning of skin representations occurs at the thalamo-cortical level while 
cortico-cortical connections serve a stereotyped competition mechanism that shapes the 
receptive fields. To further assess this hypothesis and the validity of the model, we 
reproduced in this article the exact experimental protocol of DiCarlo et al. that has been 
used to examine the structure of receptive fields in area 3b of the primary somatosensory 
cortex. Using the same analysis toolset, the model yields consistent results, having 
most of the receptive fields to contain a single region of excitation and one to several 
regions of inhibition. We further proceeded our study using a dynamic competition that 
deeply influences the formation of the receptive fields. We hypothesized this dynamic 
competition to correspond to some form of somatosensory attention that may help to 
precisely shape the receptive fields. To test this hypothesis, we designed a protocol where 
an arbitrary region of interest is delineated on the index distal finger pad and we either (1) 
instructed explicitly the model to attend to this region (simulating an attentional signal) (2) 
preferentially trained the model on this region or (3) combined the two aforementioned 
protocols simultaneously. Results tend to confirm that dynamic competition leads to 
shrunken receptive fields and its joint interaction with intensive training promotes a 
massive receptive fields migration and shrinkage. 
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INTRODUCTION 

In a previous work (Detorakis and Rougier, 2012), we proposed a 
computational model of the somatosensory cortex based on neu- 
ral field theory (Amari, 1977; Bressloff, 201 1). This model allowed 
us to investigate formation and maintenance of ordered topo- 
graphic maps in the primary somatosensory cortex during the 
critical period of development (postnatal), where representations 
are shaped, and the post-critical period, where representations are 
maintained and possibly reorganized in face of cortical or sen- 
sory lesions or dynamic changes of the environment. The main 
hypothesis of the model is that feed- forward thalamocortical con- 
nections are an adequate site of plasticity while cortico-cortical 
connections drive a competitive mechanism that is central in 
the learning process. The model relies functionally on the bal- 
ance between lateral excitation and inhibition, allowing to widen 
or sharpen the response of the model and plays a critical role 
in the shaping of the receptive fields during development. This 
modulation of the balance may originate from at least two dis- 
tinct processes at two different time scales. In the long-term. 



neurogenesis/neuronal death and synaptogenesis/synaptic degen- 
eration (Edelman, 1987) are ontogenetic factors that shape corti- 
cal connectivity during development as explained in Bressler and 
Tognoli (2006). Synaptic density spikes during the childhood fol- 
lowed by a decline during adolescence and adulthood (Feinberg 
et al, 1990). 

To further support this hypothesis, we first reproduced in this 
article the experimental protocol of DiCarlo et al. (1998) that has 
been used to characterize the structure of receptive fields (RFs) 
in area 3b of primary somatosensory cortex in three alert mon- 
keys. This protocol is based on the passive stimulation of the 
distal finger pad using a rotating drum. This allowed the authors 
to show that most RFs contain a single, central region of exci- 
tation and one or more regions of inhibition. In this work, we 
adapted this protocol to our model and validated our results 
using the same modified linear regression algorithm to char- 
acterize excitatory and inhibitory components of each RF. This 
helped us to tune the model and we found very consistent results 
using a stereotyped profile for lateral connections, resulting 
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from a fixed balance between the amount of excitation and 
inhibition. 

We further processed our analysis by considering the dynamic 
modulation of the competition following a top-down signal that 
is supposed to originate from higher order cortical areas and has 
been implemented as a gain multiplication at the level of the 
lateral intra- cortical connections. In the short-term perspective, 
such modulation allows the model to give a sharper and stronger 
response to any stimulus. In the long term perspective, the 
repeated modulation of the response has a long-lasting influence 
onto the structure of the RFs. We hypothesized such a modulation 
to represent a form of somatosensory attention (spatial atten- 
tion) because such modulation has been already proposed in the 
visual dimension as a possible mechanism for spatial attention, 
more specifically in area V4 (Salinas and Abbott, 1997; Salinas 
and Sejnowski, 2001). Indeed, attention has been mostly stud- 
ied in the visual system and can be defined as a mechanism 
that enhances the processing of interesting (understood as behav- 
iorally relevant) locations (spatial or featural) while darkening 
the rest (Posner, 1980; Treisman, 1988). The first neural cor- 
relate of that phenomenon has been discovered by Moran and 
Desimone (1985) in V4 where neurons respond preferentially 
for a given feature in their receptive fields. Since then, atten- 
tional effects have been found in each map of the ventral stream 
but also in the dorsal stream (area MT encoding for stimulus 
movement, LIP representing stimuli in a head-centered reference 
frame). Such attentional effects have also been identified in other 
modalities as well: auditory (Picton and Hillyard, 1974; Fritz et al., 
2007), motor (Norman and Shallice, 1980) and somatosensory 
to a much lesser extent (Hsiao et al, 1993). In fact, even if the 
somatosensory system has been extensively studied in monkeys 
and rats, the nature of attentional mechanisms and how they may 
affect neocortical maps of somatosensory cortices remain largely 
unknown. 

Our main hypothesis is that the modulation of a response in 
area 3b may be one of the core mechanism, even though the ori- 
gin of the modulation signal is not detailed in this article. To test 
this hypothesis, we developed a specific protocol where modula- 
tion occurs only if a presented stimulus is located within a region 
of interest (Rol) that corresponds to the attended region and we 
compared results with a protocol where the region of interest is 
specifically trained. Results tends to highlight a prominent role of 
the modulation into the shrinkage of the RFs even if only the joint 
interaction of training and attention lead to maximal effects. 

MATERIALS AND METHODS 

MODEL 

Finger pad 

We modeled a skin patch of the index distal finger pad where 
Merkel's ending complex (MEG) density is known to be the high- 
est and to convey information about touch and pressure (Pare 
et al., 2002). These receptors have been shown to have a sustained 
response to any mechanical deflection of the skin tissue. We thus 
considered a set of 256 receptors uniformly spread over the skin 
patch. When a stimulus is applied at a given position z of the skin 
patch, its mechanic property extends the pressure level to nearby 
locations (Goodwin et al, 1995). More formally, the response 



of any receptor i located at ri is given by the following equation: 

Si(z) = exp|^-^||z- rill^ (I) 

It is apparent that when a stimulus is present and its distance from 
the corresponding receptor tends to zero, the activity is the high- 
est possible. On the contrary, when there is no stimulus present, 
the activity is zero. This model assumes a very simple correlation 
between the distance of the receptor to the stimulus center and 
its level of activity. We chose such a simple model because it eases 
the mathematical analysis of the model and we are not interested 
in the full modeling of the finger pad. More accurate models can 
be found in Srinivasan (1989) (waterbed model), Dandekar et al. 
(2003) (finite elements) and in Sripati et al. (2006b) (continuum 
mechanics) but we do not think using these models would fun- 
damentally change the properties of our model (see Figure 1 for 
a comparison of the waterbed and Gaussian surface deflection 
models) since the set of 256 receptors encode a two-dimensional 
quantity that corresponds to the position of the stimulus. 

Dorsal patliway 

The dorsal column-medial lemniscus (DGML) pathway is the 
major afferent pathway for mechanosensory information and 
mediate tactile discrimination as well as proprioception (Purves 
et al., 2001). There exist several relays along this path (dorsal root 
ganglion, gracile and cuneate nuclei of caudal medulla and ventral 
posterior lateral nucleus of the thalamus) that convey informa- 
tion from first order neurons up to the somatosensory cortex. 
We modeled this complex pathway as a direct transformation of 
the MEG activity corresponding to the mean distance between 
receptors activity and the corresponding feed-forward weights. 
Gonsequently, and considering a stimulus at position z on the skin 
patch, the input /(x, z, t) received by a neuron x of SI is given by 
equation: 

1 " 

I{x,z,t) = \--T\si(z)-wUx,t)\ (2) 
n ^ 

i = 0 

where i designates a specific skin receptor and w|(x, t) is the feed- 
forward weight at time t linking receptor i to neuron x. This 
equation implies that any SI neuron receives input from all the 
skin receptors. From a neurophysiological point of view, such 
an assumption is valid to the extent that we considered only a 
small skin patch on distal finger pad. The transformation itself 
can be considered as the complement of the normalized distance 
between the set of receptors and the set of feed-forward weights. 
Such transformation is maximal (/(x, z) = 1) for a given stimulus 
z if Vi, Sf(z) = w|(x, t). This is true because Equation (1) implies 
that the maximum amplitude of a stimulus is equal to one and we 
assumed that the feed-forward weights, w^, are bound between 
0 and 1 and therefore the maximal value of /(x, z) = 1 and the 
minimum value can be /(x, z) = 0. 

Area 3b 

Area 3b of the somatosensory cortex has been modeled using 
neural field theory (Wilson and Go wan, 1973; Amari, 1977; 
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FIGURE 1 I Skin model. The finger pad skin patch is approximately of size 
25 mm^, using a receptor density of 10 mm^. It has been modeled as a planar 
surface and we considered 256 MEC's that are arranged in a regular grid over 
the whole surface with a position jitter of 5%. This results in a quasi-uniform 
distribution consistent with actual distribution of MEC as reported in Pare et al. 



(2002). (A) Schematic diagram of the hand. (B) Position and relative size of the 
skin patch. (C) Magnification of the skin patch showing MECs distribution. (D) 
Waterbed surface deflection model from Srinivasan (1989). (E) Gaussian 
surface deflection model from Detorakis and Rougier (2012). Each model 
predicts smaller deflection as a function of the distance from the load. 



Taylor, 1999) which considers the cortex as a continuous surface 
Q. Considering a stimulus z, the dynamic of the field is given by 
equation: 

decay term ' v ' feed-forward input 

lateral interaction 

(3) 

where u(x, t) is the membrane potential at position x, r 
is the membrane time constant, / is the firing rate function, w/ 
is the lateral connections function and /(x, z, t) is the output 
from the DCLM pathway as defined in previous section (see 
Figure 2). The dynamic of the field is tightly linked to the lateral 
connections function w/ that defines the behavior of the field 
(traveling waves, spiral waves, bump solutions, see Bressloff 
(20 11) for extensive review). In Detorakis and Rougier (2012), we 
defined w/ as a difference of Gaussian functions such as to obtain 
bump solutions. More precisely, we assume w/ is both isotropic 
and homogeneous (i.e., w/(x, y) = w/(|x — y|)) and defined as 

wi{x) = We{x) - Wiix) = Ke exp (y^) - Ki exp (-^) (4) 



excitatory part inhibitory part 



where (Ke, Gq) and {Ki, Gi) are constants that describe the extent 
and the strength of short-range excitation and long-range inhibi- 
tion [oi > ag). 

Learning occurs at the thalamo- cortical level using an Oja-like 
learning rule (proportional to a pre-synaptic measure multiplied 
by a post- synaptic quantity) which solves stability problems that 
is known to exist in the standard Hebbian learning rule see Oja, 
1982). It reads: 

awf (x, t) r 

,V = / (^(z) - ^/(x, t)) / we(\x - YDfiuiY, t))dy (5) 
^ V 'l^ 

pre-synaptic term ' "-^ ' 

post-synaptic term 

where y is a constant learning rate. We showed in Detorakis and 
Rougier (2012) how this learning rule, coupled with the neural 
field, allow the model to self- organize and develop topological 
representations of the skin patch. All the details are given in 
Detorakis and Rougier (2012) but briefly. Equation (3) allows 
the model to exhibit a single bump of activity (for any input) 
and the learning rule (5) exploits this bump solution to promote 
learning at position where the excitatory part of the lateral con- 
nections function is maximal. It is to be noted that because of 
the pre-synaptic term and the boundedness of receptors values 
(i.e., are bounded between 0 and 1), feed-forward weights are also 
bounded between 0 and 1. 
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Gain modulation 

As explained earlier, the shape of the bump solution of the neu- 
ral field can be controlled via lateral connections function w/. 
We have been using until now a stereotyped profile defined by 
the extent and the strength of short-range excitation (Ke, a^) and 




Cortical Layer 



Neuron 



Receptor 



Stimulus 



Input Layer 



FIGURE 2 I Schematic of the full model. Area 3b has been modeled using 
a neural field with lateral short-range excitation (wg) and long-range 
inhibition (wj). Each unit is fed with the information from all the 256 MEC 
receptors via feed-forward connections (wf). 



long-range inhibition (Ki, Oi). This profile is used for the whole 
duration of the initial training protocol and has a direct influence 
on the self- organization process. We could have used instead a 
wider/weaker or thiner/stronger profile as shown in Figure 3 but 
more importantly, we can also modify it online, provided a sig- 
nal is sent to indicate which profile is to be used for processing 
the next stimulus. This is what we refer as the attentional signal, 
originating from higher cortical areas. More precisely, we can use 
two parameters sets, (K^, K-) and {K'^ , K'^), and use the first set 
when no attentional signal is present and the second one, when 
an attentional signal is present. 

PROTOCOLS 
Initial training 

Since the model initially possesses random weights, it is firstly 
necessary to train it in order to develop topological represen- 
tations of the skin patch. We thus re-implemented the training 
protocol that has been used in Detorakis and Rougier (2012) and 
the training set is made of 50000 stimuli with random positions 
uniformly distributed over the whole skin patch. Each stimulus is 
presented once to the model and Equations (3) and (5) are eval- 
uated simultaneously until stability is achieved, i.e., there is no 
noticeable difference between u{t) and u{t -\- 8t). The model is 
reset and another stimulus is picked up until there is no more 
available stimuli. A significant difference with the original model 
is the toric implementation of skin patch and cortical model 
as well. This means that any part of a stimulus that lay out- 
side the skin patch reappears at the opposite side of the skin 
patch. The same holds true for the cortical sheet. This has been 
done to avoid any boundary effects that are known to exist in 
self- organization models. Once the training ends, the model has 
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FIGURE 3 I Gain modulation. The response of the nnodel depends 
functionally on the balance between lateral excitation (gain Kq) and 
inhibition (gain Kj), allowing to widen (A) or sharpen (C) the peak of 
activity when a stinnulus is presented. If we consider the trigger threshold 



to be the peak of nonninal response (B), the same stinnulus can either 
trigger a sharp response or not trigger any response at all, depending on 
the nnodulation. This modulation is considered in this work as a form of 
somatosensory attention. 
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developed a topological representation of the skin patch such that 
two neighbor neurons on area 3b represent two neighbor location 
on the skin patch. 

Drum protocol 

The drum protocol is a direct adaptation of the protocol that has 
been used in DiCarlo et al. (1998). Authors used a cylindrical 
drum covered with a plastic sheet (28 x 250 mm) that possesses 
raised dots pattern (with a density of 10 dots per square centime- 
ter for a total of 750 dots). The drum was mounted on a rotating 
drum stimulator and the orientation and the angular velocity of 
the drum were adjusted to produce proximal-to-distal stimulus 
movement at 40 mm/s across the skin surface. The drum com- 
pleted 100 revolutions and the total time of simulation was 14 min 
and stepped a total distance of 20 mm. We adapted the drum pro- 
tocol as a planar surface of size 250 x 30 mm and moved the skin 
patch over the full length (40 mm/s) before jumping back to the 
start and shifting up the patch by 200 /xm. The drum surface is 
made of 750 uniformly distributed dots, achieving a mean den- 
sity of 10 dots/cm^. Using a sample time step of 5 ms, the model 
has been fed with 120000 samples for a complete sweep of the 
drum surface. Activity of all neurons are recorded at once without 
centering the drum onto each individual receptive field. 

Rol protocol 

We first defined an arbitrary region of interest (Rol) on the sur- 
face of skin patch whose size is one quarter of the total skin patch 
surface (see Figure 4C, the shaded squared area in the middle of 
the skin patch). For the intensive training session, we used a set 
of 25000 stimuli such that one out of two stimuli landed into the 
Rol [1 in / 1 out ratio, (see Figure 4)]. This means that the Rol, 
was twice more stimulated compared to the rest part of the skin 
patch. We presented each stimuli once to the model until no more 
stimuli were available. Learning occurs for the whole duration 
of the protocol. For the attentional experiment, we used 25000 



uniformly spread stimuli over the whole skin patch. We presented 
each stimulus once to the model until no more stimuli were avail- 
able. If a stimulus position was within the Rol, (i.e., the center 
of the stimulus, which is the most active zone of a stimulus) we 
explicitly instructed the model to attend to this stimulus by mod- 
ifying the gain of the lateral connections (Ke and Ki) as explained 
in the gain modulation section. Learning occurs for the whole 
duration of the protocol. 

RESULTS 

CHARACTERIZATION OF THE RFs 

We first report results concerning the characterization of RF 
structures observed in area 3b following the exact protocol of 
DiCarlo et al. (1998), that was used to investigate the two- 
dimensional structure of area 3b neuronal non- classical receptive 
fields (ncRFs) in three alert monkeys (non-classical receptive 
fields are defined in Supplementary Material). As explained ear- 
lier, this protocol has been slightly adapted to meet the constraints 
of the proposed computational model architecture. Following 
the initial training protocol were cortical representations have 
been shaped (see Detorakis and Rougier, 2012), we applied the 
drum protocol for a total of 120000 samples (that can have 
multiple sites of skin patch stimulation because of the raised 
dot patterns). From these data, we applied the exact same lin- 
ear regression algorithm proposed and used by DiCarlo et al. 
(1998) for the characterization of the excitatory and inhibitory 
components of each ncRF. More precisely, for each matrix repre- 
senting an non-classical receptive field, we first convolved it with 
a Gaussian filter {fi = 0 and a = 1.7) and applied a threshold- 
ing (10% of the absolute peak value) on every value. If a value 
was below the threshold, it was set to zero. We let each pixel 
of the non- classical receptive field to have at least two of the 
four neighbors non-zero and of the same sign such that isolated 
islands of positive or negative values were not allowed if they 




mm mm mm 



FIGURE 4 I Protocol stimuli sets. (A) The training protocol set is nnade of 
50000 stinnuli distributed uniformly over the whole skin patch. At any 
moment, only one stimulus is presented to the model. (B) The drum 
protocol is based on a rotating drum made of 750 dots spread over the 
surface of the drum. The rotation of the drum makes stimuli to enter on 
the left side and exit on the right side of the skin patch, leading to 
temporal correlation between the different trials. At any moment, one to 



several stimuli can simultaneously stimulate the skin patch. (C) The Rol 
protocol, in the case of intensive training, is made of two sets of equal 
size (12500) for a total of 25000 stimuli. One set is made of stimuli 
exclusively located in the center of the skin patch and the other set is 
made of stimuli located outside this central region. This results in a higher 
(twofold) stimulus density in the central region. At any moment, only one 
stimulus is presented to the model. 
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had a total area less than 0.7 mm^. Each time we computed a 
ncRF, we also computed the signal-to-noise ratio (SNR) as well 
as the noise index, in order to constraint them to low values (see 
Supplementary Material). After this preprocessing stage, we mea- 
sured the respective size of excitatory (positive) and inhibitory 
(negative) areas. The minimum and maximum values of excita- 
tory ncRFs were 9.12, 25.92 mm^, respectively for a mean size 
of 14.142. jYiQ minimum and maximum values of inhibitory 
ncRFs were 5.92, 26.56 mm^, respectively for a mean size of 
14.4^. Figure 5 shows the bivariate plot of excitatory vs. inhibitory 
area (similar results have been found by DiCarlo et al., 1998). 
Furthermore, a k-means classification of the ncRFs was per- 
formed on the ncRFs in order to compare the number of ncRFs 
classes from the model with the number of classes in DiCarlo et al. 
(1998). The k-means classification separated 16 different classes 
according to the topology of the excitatory and inhibitory areas 
(homogeneity = 0.39, completeness = 1.0, V-measure = 0.56). 
We found non-classical receptive fields whose excitatory area was 



surrounded by the inhibitory one as well as non-classical recep- 
tive fields whose excitatory area was facing the inhibitory area (see 
Figure 5). It is to be noted that Figure 5 shows a remarkable sim- 
ilarity with physiological results of DiCarlo et al. where most of 
the ncRFs are centered around a central point of 15mm2 (exci- 
tatory) / 15 mm^ (inhibitory). The spread is larger in the case of 
DiCarlo but this was expected since we used a toric stimulation. 

TRAINING THE Rol 

During the specific training of the Rol, we considered a set of 
25000 stimuli, half of them being located in the Rol. We will later 
refer to this as the intensive protocol. At the end of the protocol, 
we measured the location and the size of the classical receptive 
fields (cRF) or simply receptive fields (RF) (see Supplementary 
Material for details) and compared them to the control setup, 
that corresponds to the end of the nominal training period (or 
the start of this protocol). Figure 6B reveals a strong migration of 
most RF toward the Rol with an overall final density being higher 
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FIGURE 5 I Characterization of the ncRFs. From the experimental drum 
protocol of DiCarlo et al. (1998), we recorded 120000 responses for each 
of the 1024 neurons of the model and we subsequently applied the 
same analysis with DiCarlo et al. (1998) to obtain the respective ncRF. 
The scatter plot on the right displays the balance between excitatory and 
inhibitory components of each ncRF. Excitatory area was measured as 
the total positive area in the thresholded ncRF (positive ncRF regions 
with values >10% of the peak absolute ncRF value, see Materials and 
Methods and Supplementary Material). Inhibitory area was measured as 



the total negative-thresholded ncRF area (negative ncRF regions with 
absolute values <10% of the peak absolute ncRF value). The left part of 
the figure (A-P) illustrates the diversity of ncRFs and the letter 
corresponds to a point in the scatter plot. The bottom row shows the 
distributions of the sizes of ncRFs. The y-axis indicate the number of 
neurons (n = 1024) and the x-axis, from left to right displays the 
excitatory area of ncRFs, the inhibitory area, the total area (is the sum 
of the excitatory and inhibitory areas) and the ratio of excitatory area to 
inhibitory one in logarithmic scale. 
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FIGURE 6 I RF Migrations. (A) The distribution of RFs over the skin patch is 
quasi-uniformly distributed for the control. (B) Intensive training onto the Rol 
makes RFs to migrate toward the Rol leading to a higher density of RFs 
within the Rol. (C) Explicitly attending the Rol modifies only marginally the 



distribution of RFs that tend to remain quasi-uniformly distributed over the 
whole skin patch. (D) The joint effect of intensive training and attention leads 
to an even greater migration of RFs toward the Rol (compared to intensive 
training only). 



in the center of the Rol in contrast to normal case illustrated in 
Figure 6A. We also measured RFs size at the end of the proto- 
col and compared them with control. The control setup shows 
a normal distribution of sizes around a central value (2.1 mm^, 
SD = 0.42) while the intensive training setup leads to a signifi- 
cant reduction of the RFs (1.6 mm^, SD = 0.48). Overall, there 
has been a significant decrease in the mean size of RFs (see 
Figures 7B\B^^ for an isolated RF and compare with the normal 
case in Figures 7A^A^^ respectively). Such results are consistent 
with Xerri et al. (1994) that shows that intensive training over a 
skin area can cause the corresponding cortical territory expan- 
sion with a simultaneous shrink of receptive fields of neurons of 
the somatosensory cortex. 

MODULATING THE Rol 

In order to make the model to attend to the Rol, we consid- 
ered a set of 25000 stimuli, uniformly spread over the whole 
skin patch and we instructed the model to attend to a stimu- 
lus if this was within the Rol, i.e., using different gains for the 



lateral connections. The major difference compared to the inten- 
sive training experiment is the non-migration of the RFs toward 
the center of the Rol as shown in Figure 6C. The distribution 
remains actually quasi- uniform and the Rol does not benefit from 
significant higher density. However, the sizes of the RFs have 
shrunk by 33%, leading to a mean size of 1.4 mm^ {SD = 0.37). 
In addition. Figures 10! show the histogram of shrinkage and 
the shrinkage of an individual RF, respectively. This demonstrates 
that migration and shrinkage of RFs are actually two distinct 
processes that can be (partly) separated. 

JOINT EFFECT OF TRAINING AND MODULATION 

For studying the joint effect of training and modulation, we 
mixed the two Rol protocols and considered both a non-uniform 
set of 25000 stimuli, half of them being located in the Rol and 
we instructed the model to attend to a stimulus if it was located 
in the Rol. The final density of RFs shown in Figure 6D reveals 
a massive migration of the RFs toward the Rol with a simultane- 
ous shrinkage in their sizes compared to the control conditions 
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FIGURE 7 I RF Shrinkage. (A') The relative histograms of RFs sizes after 
initial training (50000 samples) follows a normal distribution. (B') After model 
training, specifically in the Rol (with a 1/1 ratio) using 25000 extra samples, 
the mean RF size has been reduced by 25% compared to the nominal mean 
size. iO By sharpening the model response when a stimulus is presented 
within the Rol (25000 samples), the mean RF size has been reduced by 33% 



compared to the nominal mean size. (D') The joint effect of training and 
modulation (25000 samples) leads to a dramatic shift in relative size of RF, 
with a mean size being half of the nominal mean size. (A''-D'') Receptive 
field of a single cell recorded at the end of each of the aforementioned 
experiments. The receptive field size in the attentional/intensive condition 
(0.007 mm^) has shrunk to one third of the control size (0.024 mm^). 



(0.71 mm^, SD = 0.04). These results point out that the com- 
bined effects of intensive training and modulation actually sum 
up, leading to both a massive migration and a dramatic shrinkage 
of RFs, down to half the nominal size (see Figures 7D^D^0• 

DISCUSSION 

Using the model presented in Detorakis and Rougier (2012), we 
first validated it using the protocol and neurophysiological data 
from DiCarlo et al. (1998). We adapted the protocol to meet com- 
putational constraints and relevant recorded data. Results clearly 
indicate that the model is able to capture the main aspects of 
the original data recorded on three alert monkeys with most 
of ncRFs to contain a single region of excitation and one or 
more regions of inhibition located on one, two, three, or all 
four sides of the excitatory center. This is the first, to the best 
of our knowledge, computational model of area 3b that is able 



to replicate real neurophysiological data with such accuracy even 
though we used a very simple model for the distal finger pad and 
the dorsal column-medial lemniscus, as well. This tends to con- 
firm that the thalamo- cortical feed- forward connections are an 
adequate site of plasticity while cortico- cortical connections drive 
the competition mechanism. Furthermore, even if the present 
study has been circumscribed to the spatial characteristics of the 
receptive fields, Sripati et al. (2006a) have shown spatio-temporal 
receptive fields (STRF) in area 3b tend to have early excitatory 
region followed by in-field (replacing) greater inhibition. Authors 
conclude that such greater inhibition observed in cortical STRFs 
points to the existence of underlying intracortical mechanisms that 
is very consistent with our own hypothesis. To go further in this 
direction, we would need to consider finite transmission speed in 
cortico -cortical connections instead of instantenous connection 
(Hutt and Rougier, 2010). 
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We have also shown how this competition mechanism can 
be expHcitly modulated by the modification of the gain at both 
the excitatory and inhibitory lateral connection levels. Such 
instructed modulation leads to receptive fields shrinkage in the 
region of interest while keeping intact the overall organization, 
with no noticeable migration of RFs. We identified such modula- 
tion as a form of spatial attention that is believed to be deployed 
selectively on this or that part of the body. Interestingly enough, 
these effects are known to occur in the visual system and a num- 
ber of recent studies have identified such effects in area MT 
(Womelsdorf et al, 2008; Anton-Erxleben et al, 2009). More pre- 
cisely, authors have shown how attention inside the cRF shrinks 
it, whereas directing attention next to the cRF expands it. Authors 
uses in their modeling work a bell shaped attentional signal while 
we have been using a constant attentional signal, modifying the 
gain for the whole population at once. 

However, in the literature, the evidence for the effects of such 
spatial attention on SI are still contradictory. Hsiao and Vega- 
Bermudez (2001 ) has shown that attention is engaged in the mod- 
ification of RFs in both primary and secondary somatosensory 
cortices and Braun et al. (2002) have confirmed such engage- 
ment of attention in the primary somatosensory cortex using 
neuroimaging techniques. However, Godde et al. (2000) claimed 
that attention is not critical in enhancing performance during a 
discrimination task even though, consecutive training or pair- 
ing stimulation (leading to co- activation) can affect the RFs. In 
the case of intensive training of the Rol, our model tends to 
suggest a large expansion of the cortical territory with a simulta- 
neous shrinkage of the receptive fields as well as strong migration 
of their centers toward the Rol. These findings are still con- 
tradictory with psychophysical and neurophysiological studies 
such as Recanzone et al. (1992); Godde et al. (2000); Pilz et al. 
(2004), where authors noticed that the cortical representations 
undergo an expansion but at the same time RFs undergo a similar 
expansion. However, other neurophysiological and neuroimaging 
studies have shown that when cortical representations expand, 
RFs sizes seem to decrease (Xerri et al, 1994; Elbert et al, 1995). 
These latter results are also consistent with early findings of Sur 
et al. (1980). They found, from neurophysiological recordings and 
mappings that the magnification factor of cortical representations 
is related to the size of RFs. More precisely, the magnification fac- 
tor is proportional to the size of RFs (the smaller the RFs the 
larger the cortical representation). Our findings tend to confirm 
that cortical representations in the case of intensive stimulation 
increase their relative size with a simultaneous RFs shrinkage. 
These findings indicate that there are two distinct processes at 
work, namely modulation and training, that are believed to be 
present simultaneously in most cases, while there may exist a few 
cases where only one process is active. This may reconcile the 
aforementioned contradictory results. To confirm these findings, 
it would thus be necessary to setup new experiments where mod- 
ulation and training needs to be carefully dissociated. This can be 
done, for example, by precisely controlling the amount of train- 
ing received by a subject and by distracting the subject such as 
drifting attentional process away from the primary task. 

Even if our model suggests a hypothesis on how somatosen- 
sory spatial attention may modify the processing of stimulus 



and promote reshaping of RFs, nothing has been said so far 
about the exact nature, the origin and the selectivity of such 
attentional signal. Sarter et al. (2005) have proposed a possible 
circuitry involving the basal forebrain corticopetal cholinergic 
system since it has been observed in several studies (Donoghue 
and Carroll, 1987; Jimenez- Capdeville et al., 1997) that the loss of 
cortical cholinergic system directly impacts attentional functions. 
Furthermore, Juliano et al. (1991) have shown that the choliner- 
gic depletion prevents expansion of somatosensory topographic 
maps, suggesting that cholinergic neurotransmitters are critical 
in the structure of cortical representations. Similarly, Rasmusson 
and Dykes (1988); Tremblay et al. (1990a,b) proposed that a 
cholinergic signal is responsible for the gain modulation of neu- 
ronal populations and that the co- activation of basal forebrain 
and the somatosensory cortex by cutaneous stimulation lead 
to enhanced cortical activity. Overall, such a signal may origi- 
nate from a complex network involving the insular cortex, the 
dorsolateral prefrontal cortex, the posterior parietal cortex, the 
ventromedial prefrontal cortex, the posterior cingulate cortex and 
the anterior cingulate cortex as proposed by Menon and Uddin 
(2010). The main point is that the insular cortex acts as a switch 
between two different prefrontal networks leading to an atten- 
tional effect through saliency occurring in the anterior insular 
cortex. We can thus speculate that such a cholinergic signal may 
affect the gain of intra- cortical lateral connections and the explicit 
signal that has been used during the attended Rol protocol may 
originate from a frontal decision. 

Finally, even though we hardly notice it in our everyday life, 
somatosensory attention plays a critical role in our perception of 
the outer world. For example, the contact of clothes on the skin 
can be largely unattended even though all body receptors are acti- 
vated at once. This results from habituation and yet, it is still pos- 
sible to concentrate on a specific part of the body to actually expe- 
rience the contact. Such spatial selectivity is very similar to the 
concept of the spotlight of attention proposed by Posner (1980) in 
the eighties for the visual perception. At that time, authors were 
hypothesizing for the existence of a dedicated control mechanism 
even though this view was later challenged by the premotor theory 
of attention proposed by Rizzolatti and Craighero (1988). This 
later theory postulates instead that there is no need for two dif- 
ferent mechanisms (attention and action) and has received sup- 
port from several electrophysiological and brain imaging studies. 
However, how this theory can be adapted to somatosensory atten- 
tion remains unclear. Our model cannot answer the question on 
the selectivity since we only used a broad and constant modula- 
tion of the model. This choice has been made because we consider 
a small part of SI cortex where exactly one bump of activity can 
exist anytime. If we were to consider a larger part of SI, where for 
example several digit representations would co- exist, we would 
need a selective attentional signal to be able to direct gain mod- 
ulations to the relevant population involved in the representation 
of the Rol. This is quite a complex problem, since this would 
involve not only a sensory representation of several digits (sensory 
homunculus), but also a motor representation (motor homuncu- 
lus) and visual information as well. This is far beyond the scope 
of the present work but we think this might allow for a better 
understanding of somatosensory attention. 
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SUPPLEMENTARY MATERIAL 

The Supplementary Material for this article can be found 
online at: http://www.frontiersin.org/journal/10.3389/fncom. 
2014.00076/abstract 

REFERENCES 

Amari, S. (1977). Dynamics of pattern formation in lateral- inhibition type neural 
fields. Biol. Cybern. 27, 77-87. doi: 10.1007/BF00337259 

Anton-Erxleben, K., Stephan, V., and Treue, S. (2009). Attention reshapes center- 
surround receptive field structure in macaque cortical area mt. Cereb. Cortex 19, 
2466-2478. doi: 10.1093/cercor/bhp002 

Braun, C, Haug, M., Wiech, K., Birbaumer, N., Elbert, T., Roberts, L., et al. (2002). 
Functional organization of primary somatosensory cortex depends on the focus 
of attention. Neuroimage 17, 1451-1458. doi: 10.1006/nimg.2002.1277 

Bressler, S. L., and Tognoli, E. (2006). Operational principles of neurocogni- 
tive networks. Int. }. Psychophysiol. 60, 139-148. doi: 10.1016/j.ijpsycho.2005. 
12.008 

Bressloff, P. (2011). Spatiotemporal dynamics of continuum neural fields. /. Phys. 
AMath. r/zeor. 45:033001. doi: 10.1088/1751-8113/45/3/033001 

Dandekar, K., Raju, B. 1., and Srinivasan, M. A. (2003). 3-d finite-element models 
of human and monkey fingertips to investigate the mechanics of tactile sense. /. 
Biomech. Eng. 125, 682-691. doi: 10.1115/1.1613673 

Detorakis, G. I., and Rougier, N. (2012). A neural field model of the somatosen- 
sory cortex: formation, maintenance and reorganization of ordered topographic 
maps. PLoS ONE 7:e40257. doi: 10.1371/journal.pone.0040257 

DiCarlo, J. J., Johnson, K. O., and Hsiao, S. S. (1998). Structure of receptive fields 
in area 3b of primary somatosensory cortex in the alert monkey. /. Neurosci. 18, 
2626-2645. 

Donoghue, J., and Carroll, K. (1987). Cholinergic modulation of sensory 
responses in rat primary somatic sensory cortex. Brain Res. 408, 367-371. doi: 
10.1016/0006-8993(87)90407-0 

Edelman, G. (1987). Neural Darwinism. New York, NY: Basic Books. 

Elbert, T., Pantev, C, Wienbruch, C, Rockstroh, B., and Taub, E. (1995). Increased 
cortical representation of the fingers of the left hand in string players. Science 
270, 305-307. doi: 10.1126/science.270.5234.305 

Feinberg, 1., Thode, H. J., Chugani, H. T., and March, J. (1990). Gamma dis- 
tribution model describes maturational curves for delta wave amplitude, cor- 
tical metabolic rate and synaptic density. /. Theor. Biol. 142, 149-161. doi: 
10.1016/50022-5193(05)80218-8 

Fritz, J., Elhilali, M., David, S., and Shamma, S. (2007). Does attention play a role 
in dynamic receptive field adaptation to changing acoustic salience in al? Hear 
Res. 229, 186. doi: 10.1016/j.heares.2007.01.009 

Godde, B., Stauffenberg, B., Spengler, E, and Dinse, H. (2000). Tactile coactivation- 
induced changes in spatial discrimination performance. /. Neurosci. 20, 
1597-1604. 

Goodwin, A., Browning, A., and Wheat, H. (1995). Representation of curved sur- 
faces in responses of mechanoreceptive afferent fibers innervating the monkey 
fingerpad./. Neurosci. 15, 798-810. 

Hsiao, S., and Vega-Bermudez, F. (2001). "8 attention in the somatosensory sys- 
tem," in The Somatosensory System: Deciphering the Brains Own Body Image, ed 
R. Nelson (Boca Raton, FL: CRC), 197-217. 

Hsiao, S. S., O'shaughnessy, D., and Johnson, K. O. (1993). Effects of selec- 
tive attention on spatial form processing in monkey primary and secondary 
somatosensory cortex. /. Neurophysiol. 70, 444-447. 

Hutt, A., and Rougier, N. P. (2010). Activity spread and breathers induced by finite 
transmission speeds in two-dimensional neural fields. Phys. Rev. E 82:055701. 
doi: 10.1103/PhysRevE.82.055701 

Jimenez- Capdeville, M., Dykes, R., and Myasnikov, A. (1997). Differential control 
of cortical activity by the basal forebrain in rats: a role for both cholinergic and 
inhibitory influences. /. Comput. Neurosci. 381, 53-67. 

Juliano, S., Ma, W., and Eslin, D. (1991). Cholinergic depletion prevents expansion 
of topographic maps in somatosensory cortex. Proc. Natl. Acad. Sci. U.S.A. 88, 
780-784. doi: 10.1073/pnas.88.3.780 

Menon, V., and Uddin, L. (2010). Saliency, switching, attention and control: a 
network model of insula function. Brain Struct. Funct. 214, 655-667. doi: 
10.1007/S00429-010-0262-0 

Moran, J., and Desimone, R. (1985). Selective attention gates visual processing in 
the extrastriate cortex. Science 229, 782-784. doi: 10.1126/science.4023713 



Nordlie, E., Gewaltig, M.-O., and Plesser, H. E. (2009). Towards reproducible 
descriptions of neuronal network models. PLoS Comput. Biol. 5:el000456. doi: 
10.1371/journal.pcbi.l000456 

Norman, D. A., and Shallice, T. (1980). "Consciousness and self-regulation: 
advances in research and theory," in Attention to Action: Willed and Automatic 
Control of Behavior, eds R. J. Davidson, G. E. Schwartz, and D. Shapiro (Kluwer; 
Plenum Publ), 1-18. 

Oja, E. (1982). Simplified neuron model as a principal component analyzer. /. 
Math. Biol. 15, 267-273. doi: 10.1007/BF00275687 

Pare, M., Smith, A., and Rice, F. (2002). Distribution and terminal 
arborizations of cutaneous mechanoreceptors in the glabrous finger 
pads of the monkey. /. Comp. Neurol. 445, 347-359. doi: 10.1002/cne. 
10196 

Phillips, J., Johnson, K., and Hsiao, S. (1988). Spatial pattern representation and 
transformation in monkey somatosensory cortex. Proc. Natl. Acad. Sci. U.S.A. 
85, 1317-1321. doi: 10.1073/pnas.85.4.1317 

Picton, T., and Hillyard, S. (1974). Human auditory evoked potentials, ii: 
Effects of attention. Electroencephalogr. Clin. Neurophysiol. 36, 191-200. doi: 
10.1016/0013-4694(74)90156-4 

Pilz, K., Veit, R., Braun, C, and Godde, B. (2004). Effects of co-activation on cor- 
tical organization and discrimination performance. Neuroreport 15, 2669. doi: 
10.1097/00001756-200412030-00023 

Posner, M. 1. (1980). Orienting of attention. Q. /. Exp. Psychol. 32, 3-25. doi: 
10.1080/00335558008248231 

Purves, D., Augustine, G. J., Fitzpatrick, D., Katz, L. C, LaMantia, A.-S., McNamara, 
J. O., et al. (2001). "The major afferent pathway for mechanosensory infor- 
mation: the dorsal column-medial lemniscus system," in The Somatic Sensory 
System (Sinauer Associates), 199-202. 

Rasmusson, D., and Dykes, R. (1988). Long-term enhancement of evoked 
potentials in cat somatosensory cortex produced by co -activation of the 
basal forebrain and cutaneous receptors. Exp. Brain Res. 70, 276-286. doi: 
10.1007/BF00248353 

Recanzone, G., Merzenich, M., Jenkins, W., Grajski, K., and Dinse, H. (1992). 
Topographic reorganization of the hand representation in cortical area 3b 
owl monkeys trained in a frequency- discrimination task. /. Neurophysiol. 67, 
1031-1056. 

Rizzolatti, G., and Craighero, L. (1988). "Advances in psychological science: Vol.2. 
Biological and cognitive aspects," in Spatial Attention: Mechanisms and Theories, 
eds M. Sabourin, F. Craik, and M. Robert (New York, NY: Psychology Press), 
171-198. 

Salinas, E., and Abbott, L. F. (1997). Invariant visual responses from attentional 

gain fields. /. Neurophysiol. 77, 3267-3272. 
Salinas, E., and Sejnowski, T. J. (2001). Book review: gain modulation 

in the central nervous system: where behavior, neurophysiology, and 

computation meet. Neuroscientist 7, 430-440. doi: 10.1177/10738584010 

0700512 

Sarter, M., Hasselmo, M. E., Bruno, J. P., and Givens, B. (2005). Unraveling the 
attentional functions of cortical cholinergic inputs: interactions between signal- 
driven and cognitive modulation of signal detection. Brain Res. Brain Res. Rev. 
48, 98-111. doi: 10.1016/j.brainresrev2004.08.006 

Srinivasan, M. (1989). Surface deflection of primate fingertip under line load. /. 
Biomech. 22, 343-349. doi: 10.1016/0021-9290(89)90048-1 

Sripati, A., Yoshioka, T., Denchev, P., Hsiao, S. S., and Johnson, K. (2006a). 
Spatiotemporal receptive fields of peripheral afferents and cortical area 3b and 1 
neurons in the primate somatosensory system. /. Neurosci. 26, 2101-2114. doi: 
10.1523/JNEUROSC1.3720-05.2006 

Sripati, A. P., Bensmaia, S. J., and Johnson, K. O. (2006b). A con- 
tinuum mechanical model of mechanoreceptive afferent responses to 
indented spatial patterns. /. Neurophysiol. 95, 3852-3864. doi: 10.1152/jn. 
01240.2005 

Sur, M., Merzenich, M. M., and Kaas, J. H. (1980). Magnification, receptive-field 

area, and hypercolumn size in areas 3b and 1 of somatosensory cortex in owl 

monkeys. /. Neurophysiol. 44, 295-311. 
Taylor, J. G. (1999). Neural bubble dynamics in two dimensions: foundations. Biol. 

Cybern. 80, 5167-5174. doi: 10.1007/s004220050534 
Treisman, A. (1988). Features and objects: the bartlett memorial lecture. Q. /. Exp. 

Psychol. 40, 201-237. doi: 10.1080/02724988843000104 
Tremblay, N., Warren, R., and Dykes, R. (1990a). Electrophysiological studies 

of acetylcholine and the role of the basal forebrain in the somatosensory 



Frontiers in Computational Neuroscience 



www.frontiersin.org 



July 2014 I Volume 8 | Article 76 | 10 



Detorakis and Rougier 



A model of somatosensory attention 



cortex of the cat. i. cortical neurons excited by glutamate. /. Neurophysiol 64, 
1199-1211. 

Tremblay, N., Warren, R., and Dykes, R. (1990b). Electrophysiological studies of 
acetylcholine and the role of the basal forebrain in the somatosensory cortex 
of the cat. ii. cortical neurons excited by somatic stimuli. /. Neurophysiol 64, 
1212-1222. 

Wilson, H. R., and Cowan, J. D. (1973). A mathematical theory of the functional 
dynamics of cortical and thalamic nervous tissue. Kybernetic 13, 55-80. doi: 
10.1007/BF00288786 

Womelsdorf, T., Anton-Erxleben, K., and Treue, S. (2008). Receptive field shift and 
shrinkage in macaque middle temporal area through attentional gain modula- 
tion./. Neurosci 28, 8934-8944. doi: 10.1523/JNEUROSCI.4030-07.2008 

Xerri, C, Stern, J., and Merzenich, M. (1994). Alterations of the cortical rep- 
resentation of the rat ventrum induced by nursing behavior. /. Neurosci. 14, 
1710-1721. 



Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 06 January 2014; accepted: 29 June 2014; published online: 28 July 2014. 
Citation: Detorakis GI and Rougier NP (2014) Structure of receptive fields in a com- 
putational model of area 3b of primary sensory cortex. Front. Comput. Neurosci. 8:76. 
doi: 10.3389/fncom.20 14.00076 

This article was submitted to the journal Frontiers in Computational Neuroscience. 
Copyright © 2014 Detorakis and Rougier. This is an open-access article distributed 
under the terms of the Creative Commons Attribution License (CC BY). The use, dis- 
tribution or reproduction in other forums is permitted, provided the original author (s) 
or licensor are credited and that the original publication in this journal is cited, in 
accordance with accepted academic practice. No use, distribution or reproduction is 
permitted which does not comply with these terms. 



Frontiers in Computational Neuroscience 



www.frontiersin.org 



July 2014 I Volume 8 | Article 76 | 11 



